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CN ■ I. OVERVIEW 

A. Model of the Refractive Index Structure Function 

This paper describes modeling of variations of refractive indices n based on their structure functions T) 
^1 I?„(Ar) = C^IArr (1) 

Q . with Kolmogorov power index 

^ : 7 = 2/3, (2) 

sampled at points which are Ar apart and located in an air volume between the light source and some input pupil of 
\ a receiver. The familiar assumptions of this model are 

^ ^ • "single layer" homogeneity: the refractive index structure constant is indeed constant over the entire volume 

[>«^ \ that is probed by the electromagnetic wave, 

o> ' . 

^■f-^ . • isotropy: the structure function depends on the absolute value of the vector Ar, but not on its direction. 

In the conjugated Fourier variable / of wave numbers, the turbulent spectrum is represented by the (again isotropic) 
, Wiener spectrum $„(/) [2l|, [l^j- The expectation value of the mean appears as a constant in the Fourier kernel, and 
' the expectation value of the covariance as a cosine: 

> 

I [The imaginary part of the factor exp(27ri/ • Ar), which is induced by the shift theorem Fourier transforms and which 
represents the contributions with odd parity, vanishes by isotropy after integration.] The power spectrum is 0, 0] 



I5„(Ar) = 2 j $„(/)[! - cos(2^/ • ^v)]cPf. (3) 



1>n(/) = -^"'-'i^j^Clr'' ^ c^Clr'"'; « 0.00969315, (4) 

2r(^) 



r(/3/2) 

2 

at a power spectral index 

/3 = 7 + 3 (=11/3). (5) 
Where the variable k — 27r/ is used in lieu of /, the constant (27r)^c$ w 0.033005 takes the role of c$. 



*URL: |http: //www. strw. leldenuiilv.iil/~mathar[ Electronic address: lmathar@strw.leidenumv.nll 
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B. Circular (Projected) 2D Case 

Above, the structure function and its spectrum are defined as functions of the 3D distances Ar and a. Integration 
along two paraUel fines of sight through a layer of vertical thickness h at an air mass 1/ sin a at a projected baseline P 
for momentum number k (equal to 27r divided tfirough the wavelength of the electromagnetic wave) defines 2D phase 
structure functions 



7+1. 



h 



sma 



where 



5(7) = 



- r(-(i + 7)/2) 
r(-7/2) 



2.9144, (7 = 2/3) 



(6) 



(7) 



is a geometric constant from the integration [10, llai ISa, [23| . The information sampled along the pencils of length 
hj sin a through the layer is reduced to T)^^ a function of the two components of P in the input pupil of the receiver. 
Noll's notation ^ip(y) — 0.023rjJ'^^^i/"^^/^ of the associated 2D power spectrum of the phase differences in terms of 
the Fried parameter t-q — at the wavelength 27r/fc of the electromagnetic spectrum — [21], (21)] is equivalent as we note 
that 0.023 is the product of the ad-hoc scaling parameter [ll[ 



2[(24/5)r(6/5)]^/6 « 6.883877 



by 



r(/3/2) 



27r2+'rr(l - /3/2) 



0.0033260. 



(8) 



(9) 



The association of the constants in equations ([7]) and © is summarized in the following diagram. The role of 
5(7) switches from multiplier to divisor while transiting from real to wavenumber space. 



cl 



3D Fourier Transform 



int(!gration over h 
multiply g(7) 



2.9146*2 



2D Fourier Transform 



0.00969C2 



divide 5(7) 



0.00333C2 



The intend is to show how the phase structure function is associated with integration of a refractive index structure 
function over path lengths; the Fried parameter and some of the associated geometric constants will disappear from 
the analysis of the next section, where no such path integration is performed. 

This concludes the overview of the literature. Karhunen-Loeve (KL) modes in the projected (2D) case are not 
reviewed here; we refer to the available literature [1, [l^, E, HH] • 



C. Applications 

The aim of calculations in sectionlUis to remove the constraint of working with optical path lengths integrated along 
two thin, parallel rays through the turbulent atmosphere. Simulations of anisoplanatism, multi-conjugate adaptive 
optics, and atmospheric transfer functions for radio (as opposed to optical) receivers often integrate over extended 
fields-of-view. Essentially we factorize the problem; we provide means to fill a volume with instances of a turbulent 
refractive index, but keep aside by which weighing scheme sub-volumes accumulate the instrumental signal. 

To keep mathematics simple, this volume is given spherical shape, the shape that is symmetry-adapted to the 
isotropy of the structure function. As one extreme application of the geometry one could image the center of this 
sphere at the center of the Earth, the sphere radius set to the Earth radius plus some tenth of kilometers, and then 
evaluating refractive indices only in the thin inner hull close to the sphere surface — representing the entire Earth's 
atmosphere — computed from superposition of the KL modes. 
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II. KARHUNEN-LOEVE FUNCTIONS IN THE SPHERE 
A. Spherical Basis 

The generic form of the Karhunen-Loeve eigenvalue equation for the modes K and their eigenvalues associated 
with a covariance C sampled over a spherical region of radius R is 



/ C{\r-r'\)K{r')d^r' = B^K{r). 

Jr'<R 



(10) 

One benefit earned from the constraints of the model described in section II Al is that a radial scaling 

r^Rx (11) 
maps these cases onto a single unified model inside the unit sphere: 

j^^^^ |x - xTif (x')d^x' = ^^i^(x). (12) 

The other benefit is that the spherical symmetry of the domain of integration matches the isotropy in the covariance; 
moving on to the Fourier domain of the dimcnsionlcss cr = Rf, 

K{x') = J (faK{a-)e~'^''"' ''' (13) 

and 

C{\x-x'\)= /d3aC(a)e-2-»'^-("-"') (14) 



transforms the convolution into a product, in which the eigenvalue problem is block diagonal in the angular variables 



Rt^ 



C{a)K{a) = -^K{a). (15) 



With the same argument as for the 2D case, ie, because the exponent of C oc is too large for the Kolmogorov 
spectrum to remain integrable in the large-scale limit a —^ 0, the decomposition ()15|) is not possible for functions with 
nonzero K{a) at the origin of the cr-space: the homogeneous mode, the analog to the 2D piston mode, is left out of 
the analysis. I?„ does not convey information on the mean of the refractive index. 

B. Diagonalization in Wavenumber Space 

Subject to the usual ambiguity of any choice of basis functions, we expand the K(x) in an orthogonal Zernike basis, 
suggested in Appendix |Al serializing the eigcnfunctions for each block I of the eigenvalue problem by an index p: 

K(^\a) = ^^„,.,zi7)(^) = 4^5]^„,,,z'l^™(0.,0.)i?W(a). (16) 

n n 

The reduction of (fT^ with ^ reads 

Cc,a-'3iff(a)=A^,,4')(a); K-i^J^- (17) 
Formulation of an eigenvalue problem of matrix algebra proceeds with (|A17p . 
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TABLE L /'^,/[22''^7ri/(2n + 3)(2n' + 3)] for n,n' < 12, I even. Multiply by ^(2n + 3)(2n' + 3) and in addition by 
4.986967 to obtain Z,''^/. The table is symmetric; only the non-redundant upper right triangle is shown. 
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b4 


54 


2V 


135 


2V 


27 


2V 


385 




142^:^4 


78^24 


2156^441 


1545^^9175 


7647^^8550 




20020 


1^72 


7436^4^29 


980^^655 


9802^5^550 


5637^^4050 






209209 


7^^825 


6760^^0 


3575^(1^330 


3038^^5305 








950950 


26^7^950 


2095g7|p9 


9776|^359 










2947945 


16|^^8855 
167806100 


s^WolP 

11^^1^^78100 














1803915575 



TABLE II: 



/i'i,/[2^''^7ri/(2n + 3)(2n' +3)] for ji,n' < 13, / odd. 



1 



5 



11 



13 

-^27- 



'04^6 
3^^323 



.6J^4 



i0^5 



3091^:^8350 



1787526650 



33713J71 
7215662300 



13945^2650 



9583^^2885 
4888^^6795 
15041g5|860 



2525481805 



A core integral matrix i^, is defined to tighten the notation, 

/•CO 1 

c^J^ aV-''i?W(a)i?iV(a)daEE/«,^. (19) 
The integral over the product of Bessel Functions of the form (|A16p yields [14,, 6.574.2] [l^ p. 49] 

e - c.^ V(2n + 3)i2n' + 3)(-l)("-"')/^^;^ , ^^f)^ "'^"^^ ^''^ 

for n,n' = 1,1 + 2,1 + 4, . . .. These are the matrices to be diagonalized to establish the KL basis functions, one matrix 
for each I, 

^2-^n!n''^nJ,p — ^p^lHn',l,p- (21) 
n 

The upper right corners of the symmetric = ^n'' n shown in Table [J for even and in Table [Til for odd I. The 
individual array for a specific I is the sub-array in which all rows n < I and columns n' < I are removed. A common 
factor has been split off to write them as exact fractions. If no such scaling is chosen, tables like Table Hill result. As 
argued at the end of section III Al row n = and column n' = in Table |T] are deleted prior to the diagonalization. 

As already in the 2D case, the relative strength of values on the diagonals compared to the sub-diagonals is predictive 
of how close the KL modes resemble pure Zernike functions. An explicit table of the dominant modes is given in 
Appendix |B1 sorted with respect to decreasing ; . 



TABLE III: limited to n,n' < 11. This repeats the values of Table HIl without the scale factor. 
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0.269027 
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0.000029 


0.000003 


0.000000 


3 




0.015929 


-0.004505 


0.000493 


0.000009 


0.000001 
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0.003681 


-0.001362 


0.000180 


0.000004 
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0.001320 


-0.000558 


0.000082 


9 










0.000596 


-0.000274 
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0.000311 
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C. Real-valued Basis 



In practice, the component e™"^ of the spherical harmonics in the Zernike basis (|A2p may be replaced by switching 
in the real-valued functions Q 

^/e^cos{me), ^/e^ism{m9), (22) 

where 

'i'":? (23) 

2, m > 1 ^ ' 

is the Neumann factor. This choice maintains normalizations, and maintains the 2Z + 1 fold degeneracy of the 
Karhunen-Loeve modes. 



III. SUMMARY 



This work traced the well-known Karhunen-Loeve functions of a covariance with power spectral index 5/3 over a 
circular telescope entrance pupil back to the associated functions with power spectral index 2/3 over the sphere. The 
solution of the eigenvalue integral equation for the covariance has been formulated for expansions in a 3D Zernike 
basis. Bare Zernike functions turn out to be good approximations to the large scale modes. 



Acknowledgments 

This work is supported by the NWO VICI grant 639.043.201 to A. Quirrenbach, "Optical Interferometry: A new 
Method for Studies of Extrasolar Planets." 



APPENDIX A: 3D ZERNIKE BASIS 



1. Definition 

We define an orthogonal 3D Zernike Basis ■Z'^™'' {r, 0, (j)) over the unit sphere in strict analogy to the domain of the 
unit circle 0, 0] as a product over (i) Courant-Hilbert Jacobi polynomials G to span the radial coordinate x and (ii) 
Spherical Harmonics Y"^ to cover the azimuthal and polar angles and 6: 

Z^^^ = R(^\x)Yr{eA), 0<x<l, O<0<^, O<0<2^. (Al) 

Spherical Harmonics 



are a product of e™"^ sin'™' 9 with a polynomial in cos 6 of degree / — |m|. (There are at least three different sign/phase 
conventions in the literature, of which only one is picked by this equation. Since the manuscript does not deal 

with aspects of angular momentum coupling, the specific choice is not relevant.) The transformation to Cartesian 

coordinates (x sin cos 0, a; sin sine/), a; cos 0) has been discussed in the literature [H, H^. 

Rotation of the 3D coordinate system mixes the components of YJ™, keeping / constant d, 2.1.3][1, Q. The role 
of the azimuthal number m of the 2D case [2Tj| is therefore taken over by I to set the lowest order to the radial 
polynomials jj, [13, [2^. With this constraint, demanding orthogonality with respect to the weight function in 
< X < 1 (of the Jacobian of the spherical coordinates), and setting a smoothness-constraint at the origin x = in 
the sense that the powers of x in the polynomial are all either even or odd, the radial functions need to be proportional 
to x''Ga{q,q,x^), where 



(A3) 
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FIG. 1: Radial part Wn'{x) of the 3D Zernike basis defined in ^M^, Z = or I. 



r{q + a) 



(-1) 



r(a + l) r(g + a + fc) 



r(g + 2a)' ^r(a-fc + l)A:! r(g + A:) 



r2(g + a) 



T{q + 2a)Tiq) 

are Courant-Hilbert Jacobi polynomials [1, 22.3.3] with parameters 

a={n-l)/2, q = l + 3/2 



(A4) 

y) (A5) 
(A6) 



The characteristic difference with the 2D case, see [20| and references therein, is that the parameter q has increased 
by 1/2 to cope with the higher power of the Jacobian of the Cartesian-to-spherical (as opposed to Cartesian-to- 
circular) transformation of coordinates. [This boost of the parameter remains visible in the matrix elements (I20p 
and remains visible in the Bessel Function indices (|A16[) which are translated from integer to half- integer values.] 
Ortho-normalization based on their orthogonality 1, 22.2.2] yields 



m. ^ , (^ + l)!(2n + 2)! 



2n-i(^n + I + 2)1 {^)l{n + ly: 

^/2n + 3 .y Ga{q, q, x^), 

r(a -I- q)T{a + 1) 



x^Rll\x)Rll}{x)dx = S, 



'G(„_,)/2a-l--,; + -,a.') 



(A7) 
(A8) 

(A9) 



Some examples of these polynomials are given in Table IIVI for comparison with the 2D case [2l|, Table I] . A subset is 
also visualized in Figure [T] for comparison with [13, Fig. 1] . The nodal counts are as expected for any of the classical 
orthogonal polynomials. The overall choice of signs, ie, Rn\l) > 0, is just inherited from the Jacobi polynomials, 
and carries no physical significance. 



2. Fourier Transform 



The Rayleigh expansion of the plane wave leads to the Fourier transform in terms of spherical Bessel Functions j 
P, sect. 10], 

(AlO) 



(All) 



= J d^xe^^'''-^Rll\x)Yr{0,cb)^4WYr{0a,<l>.) J x^dxji{27rax)Rll\x). 
Repeated use of the partial integration P, 11.3.3] 

x''+^C^-i{x)dx = x^'+^C^ix) + [ x''C^{x)dx 
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TABLE IV: The radial polynomials Rii\x) for Z < 4 and n < 10. 



J. n Rn^x) 

73 

2 i^/7(-3 + 5a;2) 

4 |7TT(l5 + 63a;^-70x2) 

6 ^715 (-35 + 429a;® - 693a;'' + 315 2;^) 

8 jig ^/T9 (315 + 12155 a;* - 25740 a;« + 18018 x"- 4620x2) 

10 jig 723 (-693 + 88179 a;'°- 230945 a;** + 218790 a;^- 90090 a;* + 15015 a;2) 

1 1 a;^/5 

1 3 |a;(-5 + 7a;2) 

1 5 i\/T3a;(35 + 99a;*-126a;2) 

1 7 ^7T7x (-105 + 715 a;« - 1287 a;'' + 693 a;2) 

1 9 ^^(1155 + 20995 a;**- 48620 a;^ + 38610 a;''- 12012 a;2) 

2 2 x^V7 

2 4 iTTTa;^ (-7 + 9a;2) 

2 6 |/r5a;2 (63 + 143a;' - 198a;2) 

2 8 -L7T9a;2 (-231 + 1105 a:^ - 2145a;* + 1287 a;^) 

2 10 (3003 + 33915 - 83980 x" + 72930 - 25740 x^) 

3 3 3x^ 

3 5 i^/T3x^ (-9 + 11x2) 

3 7 |7T7x^ (99 + 195 x"- 286x2) 

3 9 i72Tx=' (-429 + 1615 x'* -3315 x'' + 2145x2) 

4 4 x-'TTT 

4 6 iyT5x' (-11 + 13x2) 

4 8 i\/T9x'' (143 + 255X''- 390x2) 

4 10 i^x* (-715 + 2261 X® -4845 x'' + 3315x2) 



for any ordinary Bessel Functions C yields 

f\'dxM2naa:)J^'^ = ±{-1)' ~ ^''^ " ^ + ^^ -.+.+1(2^) = ^ ^ (2.a). 

Jo Sttct ^ (tto-)'' 2t:(7 ^ (ttct)*^ 

(A12) 

This is the building block to transform the radial polynomials of interest, ie, those constructablc from non-negative 
integers s. Pochhammer's symbol 1, 6.1.22] 

(z)„, = z(z + l)(z + 2)---(z + m-l) = r(z + m)/r(z); (z)o = 1, (A13) 

has been used to condense the notation for the falling factorials. Summation over a terms of this type yields 

^ x^3i{2nax)x'G^{q,q,x')dx ^ ^^^^ p^^^'^ (2^a)r(a + g + fc) Q . (A14) 
Gathering also the constant of normalization in (jA7p concludes 

i?«(a) ^ j\'ji{27Tax)R(^\x)dx = ^^^^ E ("i) ^ ^^^'^^^^ + (fc) (^^^^ 

= (-l)"V2^^^^ii^^, (A16) 

271(7 

orthogonal with 

(A17) 

A graphical representation of the cases Z < 3, a < 3 is in Figure [21 The leading polynomial order is Rn\(T) cx cr" as 
cr^ 0. P, 10.1.2]. 
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FIG. 2: The simplest cases of R)i'{a), < Z < 3, q < 3, as given in (KTE^ . 



In summary, (jAlOp has been factorized as 
where Parseval's theorem maintains orthogonality in wavcnumber space, 



(A18) 



(A19) 



APPENDIX B: MODE LIST 

The dominant modes Kp'' of ((T6| are Usted below, separated by empty lines, each representing 21 + 1 modes labeled 
by m. The notation consists of the value ; followed by a colon, followed by I, followed by another colon and the 

Zernike expansion. The Zernike expansion is a sum over floating point numbers Knj.p multiplied by the latter 

written as Z(n, 1) , optionally wrapped into a new line. The squares of K„^i_p add to unity in each line. The phase/sign 
convention for the eigenvectors is that Kn,i^p with maximum absolute value is positive in each line. 
Small-terms contributions with |K„,i,p| < 10~^ are not printed. 

0.2723742 : 1: -fO . 9935210*Z(1 , 1) -0 . 1132657*Z(3, 1) -fO . 0093216*Z(5 , 1) -0 . 0001473*Z(7, 1) -f0.0000128*Z(9, 1) 
-f0.0000014*Z(ll,l) 



0.0495069 : 2: -f0.9722349*Z(2,2) -0 . 2319940*Z(4,2) -fO . 0305878*Z(6 , 2) -0.0015909*Z(8,2) -fO . 0000930*Z(10 ,2) 
+0 . 0000023*Z (12,2) +0 . 0000011*Z (14 ,2) 

0.0495069 : 0: -f0.9722349*Z(2,0) -0 . 2319940*Z(4,0) -f0.0305878*Z(6,0) -0.0015909*Z(8,0) -f0.0000930*Z(10,0) 
+0 . 0000023*Z (12,0) +0. 0000011*Z (14,0) 
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0.0174576 : 3: +0.9473332*2(3 
-0.0000056*Z(13,3) +0.0000023*Z 

0.0143736 : 1: +0 . 1081177*Z(1 
+0.0005484*Z(11,1) -0.0000162*Z 

0.0080675 : 4: +0.9224733*Z(4 
-0.0000345*Z(14,4) +0.0000046*Z 

0.0059302 : 2: +0 . 2144266*Z(2 
+0.0019187*Z(12,2) -0.0001183*Z 

0.0059302 : 0: +0 . 2144266*Z(2 
+0.0019187*Z(12,0) -0.0001183*Z 

0.0043466 : 5: +0.8987736*Z(5 
-0.0000988*Z(15,5) +0.0000094*Z 

0.0030561 : 3: +0.2833231*2(3 
+0.0044522*Z(13,3) -0.0003824*Z 

0.0028433 : 1: +0 . 0312484*Z (1 
-0.0412113*Z(11,1) +0.0054932*Z 

0.0025891 : 6: +0 . 8765467*Z(6 
-0.0002124*Z(16,6) +0.0000189*Z 

0.0017869 : 4: +0 . 3306586*Z (4 
+0.0082606*Z(14,4) -0.0008873*Z 

0.0016562 : 7: +0.8558181*Z(7 
-0.0003881*Z(17,7) +0.0000358*Z 

0.0015654 : 2: -0.0803297*2(2 
+0.0731123*2(12,2) -0.0122479*2 
-0.0000014*2(22,2) 

0.0015654 : 0: -0.0803297*2(2 
+0.0731123*2(12,0) -0.0122479*2 
-0.0000014*2(22,0) 

0.0011353 : 5: -0.3644454*2(5 
-0.0133381*2(15,5) +0.0016997*2 

0.0011176 : 8: +0.8365082*2(8 
+0.0057793*2(16,8) -0.0006360*2 

0.0009624 : 3: -0.1234664*2(3 
+0.1096522*2(13,3) -0.0219805*2 
-0.0000031*2(23,3) 

0.0009269 : 1: -0.0132054*2(1 
-0.3857727*2(11,1) +0.1191289*2 
+0.0000351*2(21,1) -0.0000037*2 

0.0007862 : 9: +0.8185038*2(9 
+0.0077854*2(17,9) -0.0009640*2 

0.0007655 : 6: -0.3891945*2(6 
-0.0196069*2(16,6) +0.0028703*2 

0.0006365 : 4: -0.1595289*2(4 
+0.1485535*2(14,4) -0.0344854*2 
-0.0000076*2(24,4) 

0.0005909 : 2: -0.0387797*2(2 
-0.4675801*2(12,2) +0.1742073*2 
+0.0001105*2(22,2) -0.0000110*2 



3) -0. 
15,3) 

1) +0. 
13,1) 

4) -0. 
16,4) 

2) +0. 
14,2) 

0) +0. 
14,0) 

5) -0. 
17,5) 

3) +0. 
15,3) 

1) +0. 
13,1) 

6) -0. 
18,6) 

4) +0, 
16,4) 

7) -0. 
19,7) 

2) -0. 
14,2) 



0) -0. 
14,0) 



5) -0. 
17,5) 

8) -0, 
18,8) 

3) -0. 
15,3) 



1) -0. 

13,1) 

23,1) 

9) -0. 
19,9) 

6) -0. 
18,6) 

4) -0. 
16,4) 



2) -0, 

14,2) 

24,2) 



3153328*2(5,3) +0.0557030*2(7,3) -0.0046947*2(9,3) +0.0003314*2(11,3) 



,9173415*2(3,1) -0.3757763*2(5,1) +0.0743888*2(7,1) -0.0072975*2(9,1) 
+0.0000031*2(15,1) 

,3771734*2(6,4) +0.0818302*2(8,4) -0.0092993*2(10,4) +0.0008094*2(12,4) 



8327378*2(4,2) -0.4927135*2(6,2) +0.1320534*2(8,2) -0.0189879*2(10,2) 
+0.0000100*2(16,2) 

8327378*2(4,0) -0.4927135*2(6,0) +0.1320534*2(8,0) -0.0189879*2(10,0) 
+0.0000100*2(16,0) 

4247285*2(7,5) +0.1076118*2(9,5) -0.0151323*2(11,5) +0.0015784*2(13,5) 



7478488*2(5,3) -0.5687499*2(7,3) +0.1890304*2(9,3) -0.0349457*2(11,3) 
+0.0000310*2(17,3) 

3323629*2(3,1) +0.6982757*2(5,1) -0.5956983*2(7,1) +0.2106757*2(9,1) 

-0 . 0004989*2(15 ,1) +0 . 0000407*2(17 , 1) 

4622201*2(8,6) +0.1324035*2(10,6) -0.0219219*2(12,6) +0.0026626*2(14,6) 



,6685227*2(6,4) -0.6180700*2(8,4) +0.2423796*2(10,4) -0.0539894*2(12,4) 
+0.0000803*2(18,4) -0.0000033*2(20,4) +0.0000010*2(22,4) 

,4923231*2(9,7) +0.1559183*2(11,7) -0.0294312*2(13,7) +0.0040659*2(15,7) 



4131467*2(4,2) -0.5450977*2(6,2) +0.6585434*2(8,2) -0.2941911*2(10,2) 
+0.0014490*2(16,2) -0.0001385*2(18,2) +0.0000078*2(20,2) 



4131467*2(4,0) -0.5450977*2(6,0) +0.6585434*2(8,0) -0.2941911*2(10,0) 
+0.0014490*2(16,0) -0.0001385*2(18,0) +0.0000078*2(20,0) 



,5960909*2(7,5) +0.6491312*2(9,5) -0.2909694*2(11,5) +0.0750988*2(13,5) 
-0.0001754*2(19,5) +0.0000115*2(21,5) -0.0000017*2(23,5) 

,5168302*2(10,8) +0.1780536*2(12,8) -0.0374644*2(14,8) 
+0.0000629*2(20,8) -0.0000023*2(22,8) +0.0000010*2(24,8) 

,4534598*2(5,3) -0.4075425*2(7,3) +0.6822917*2(9,3) -0.3674242*2(11,3) 
+0.0031596*2(17,3) -0.0003567*2(19,3) +0.0000288*2(21,3) 



1534280*2(3,1) -0.4633881*2(5,1) -0.3595625*2(7,1) +0.6845747*2(9,1) 
-0.0245691*2(15,1) +0.0036299*2(17,1) -0.0004190*2(19,1) 



5369951*2(11,9) +0.1988032*2(13,9) -0.0458634*2(15,9) 
+0.0001032*2(21,9) -0.0000060*2(23,9) +0.0000014*2(25,9) 

5304668*2(8,6) +0.6673369*2(10,6) -0.3345193*2(12,6) +0.0974584*2(14,6) 
-0.0003351*2(20,6) +0.0000280*2(22,6) -0.0000032*2(24,6) 

4702764*2(6,4) -0.2875787*2(8,4) +0.6806273*2(10,4) -0.4291202*2(12,4) 
+0.0057898*2(18,4) -0.0007570*2(20,4) +0.0000753*2(22,4) 



2229036*2(4,2) -0.4693503*2(6,2) -0.1796996*2(8,2) +0.6673417*2(10,2) 
-0.0429519*2(16,2) +0.0076259*2(18,2) -0.0010468*2(20,2) 



0.0005909 : 0: -0.0387797*2(2,0) -0.2229036*2(4,0) -0.4693503*2(6,0) -0.1796996*2(8,0) +0.6673417*2(10,0) 
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-0.4675801*2(12,0) +0. 1742073*Z(14,0) -0.0429519*2(16,0) +0.0076259*2(18,0) -0.0010468*2(20,0) 
+0 . 0001105*2 (22 , 0) -0 . 00001 10*2 (24 , 0) 

0.0005719 : 10: +0.8016857*2(10,10) -0.5537242*2(12,10) +0.2182106*2(14,10) -0.0545025*2(16,10) 
+0.0100622*2(18,10) -0.0013779*2(20,10) +0.0001598*2(22,10) -0.0000118*2(24,10) +0.0000020*2(26,10) 

0.0005398 : 7: -0.4076451*2(7,7) -0.4711126*2(9,7) +0.6763013*2(11,7) -0.3731525*2(13,7) +0.1204400*2(15,7) 
-0.0269500*2(17,7) +0.0044332*2(19,7) -0.0005783*2(21,7) +0.0000570*2(23,7) -0.0000061*2(25,7) 

0.0004436 : 5: -0.1893914*2(5,5) -0.4728633*2(7,5) -0.1841723*2(9,5) +0.6623861*2(11,5) -0.4795241*2(13,5) 
+0.1880973*2(15,5) -0.0494097*2(17,5) +0.0094393*2(19,5) -0.0014015*2(21,5) +0.0001624*2(23,5) 
-0.0000171*2(25,5) 

0.0004276 : 11: +0.7859404*2(11,11) -0.5676904*2(13,11) +0.2363442*2(15,11) -0.0632823*2(17,11) 
+0.0125856*2(19,11) -0.0018812*2(21,11) +0.0002354*2(23,11) -0.0000206*2(25,11) +0.0000029*2(27,11) 

0.0004024 : 3: -0.0651981*2(3,3) -0.2733637*2(5,3) -0.4451931*2(7,3) -0.0326736*2(9,3) +0.6203023*2(11,3) 
-0.5285997*2(13,3) +0.2303104*2(15,3) -0.0658223*2(17,3) +0.0135712*2(19,3) -0.0021567*2(21,3) 
+0.0002690*2(23,3) -0.0000292*2(25,3) +0.0000019*2(27,3) 

0.0003943 : 8: -0.4215561*2(8,8) -0.4173790*2(10,8) +0.6785254*2(12,8) -0.4071763*2(14,8) 
+0.1435708*2(16,8) -0.0352331*2(18,8) +0.0064074*2(20,8) -0.0009229*2(22,8) +0.0001035*2(24,8) 

-0.0000112*2(26,8) 

0.0003932 : 1: -0.0068258*2(1,1) -0.0831481*2(3,1) -0.2897187*2(5,1) -0.4325853*2(7,1) +0.0025750*2(9,1) 
+0.6065149*2(11,1) -0.5392619*2(13,1) +0.2411059*2(15,1) -0.0703318*2(17,1) +0.0147685*2(19,1) 
-0.0023858*2(21,1) +0.0003027*2(23,1) -0.0000331*2(25,1) +0.0000022*2(27,1) 

0.0003271 : 12: +0.7711645*2(12,12) -0.5794048*2(14,12) +0.2532829*2(16,12) -0.0721245*2(18,12) 
+0.0153311*2(20,12) -0.0024758*2(22,12) +0.0003325*2(24,12) -0.0000329*2(26,12) +0.0000043*2(28,12) 

0.0003217 : 6: -0.2141615*2(6,6) -0.4666186*2(8,6) -0.0954608*2(10,6) +0.6334239*2(12,6) -0.5195267*2(14,6) 
+0.2270429*2(16,6) -0.0663391*2(18,6) +0.0141531*2(20,6) -0.0023475*2(22,6) +0.0003081*2(24,6) 
-0 . 0000351*2 (26 ,6) +0 . 0000025*2 (28 , 6) 

0.0002963 : 9: -0.4321073*2(9,9) -0.3686297*2(11,9) +0.6757850*2(13,9) -0.4369714*2(15,9) 
+0.1665011*2(17,9) -0.0443179*2(19,9) +0.0087990*2(21,9) -0.0013848*2(23,9) +0.0001726*2(25,9) 
-0 . 0000196*2 (27 ,9) +0 . 0000011*2 (29 , 9) 

0.0002873 : 4: +0.0899412*2(4,4) +0.3087998*2(6,4) +0.4055501*2(8,4) -0.0847063*2(10,4) -0.5563662*2(12,4) 
+0.5704152*2(14,4) -0.2846000*2(16,4) +0.0922604*2(18,4) -0.0215806*2(20,4) +0.0038876*2(22,4) 
-0.0005548*2(24,4) +0.0000668*2(26,4) -0.0000059*2(28,4) 

0.0002739 : 2: +0.0217345*2(2,2) +0.1326368*2(4,2) +0.3344713*2(6,2) +0.3619721*2(8,2) -0.1549178*2(10,2) 
-0.5119804*2(12,2) +0.5869219*2(14,2) -0.3108572*2(16,2) +0.1054741*2(18,2) -0.0256715*2(20,2) 
+0.0047925*2(22,2) -0.0007084*2(24,2) +0.0000877*2(26,2) -0.0000082*2(28,2) +0.0000011*2(30,2) 

0.0002739 : 0: +0.0217345*2(2,0) +0.1326368*2(4,0) +0.3344713*2(6,0) +0.3619721*2(8,0) -0.1549178*2(10,0) 
-0.5119804*2(12,0) +0.5869219*2(14,0) -0.3108572*2(16,0) +0.1054741*2(18,0) -0.0256715*2(20,0) 
+0.0047925*2(22,0) -0.0007084*2(24,0) +0.0000877*2(26,0) -0.0000082*2(28,0) +0.0000011*2(30,0) 

0.0002552 : 13: +0.7572651*2(13,13) -0.5892625*2(15,13) +0.2691084*2(17,13) -0.0809680*2(19,13) 
+0.0182744*2(21,13) -0.0031620*2(23,13) +0.0004537*2(25,13) -0.0000497*2(27,13) +0.0000062*2(29,13) 

0.0002407 : 7: -0.2348116*2(7,7) -0.4548404*2(9,7) -0.0194882*2(11,7) +0.5976955*2(13,7) -0.5502531*2(15,7) 
+0.2645281*2(17,7) -0.0848500*2(19,7) +0.0199308*2(21,7) -0.0036443*2(23,7) +0.0005316*2(25,7) 
-0 . 0000659*2 (27 ,7) +0 . 0000060*2 (29 , 7) 

0.0002280 : 10: -0.4401180*2(10,10) -0.3242846*2(12,10) +0.6693673*2(14,10) -0.4629349*2(16,10) 
+0.1889767*2(18,10) -0.0540704*2(20,10) +0.0116042*2(22,10) -0.0019775*2(24,10) +0.0002699*2(26,10) 
-0 . 0000326*2 (28 , 10) +0 . 0000025*2(30 , 10) 

0.0002127 : 5: +0.1122926*2(5,5) +0.3328253*2(7,5) +0.3585602*2(9,5) -0.1767793*2(11,5) -0.4838664*2(13,5) 
+0.5954557*2(15,5) -0.3351895*2(17,5) +0.1212944*2(19,5) -0.0316439*2(21,5) +0.0063570*2(23,5) 
-0.0010169*2(25,5) +0.0001357*2(27,5) -0.0000143*2(29,5) +0.0000018*2(31,5) 

0.0002025 : 14: +0.7441600*2(14,14) -0.5975741*2(16,14) +0.2839011*2(18,14) -0.0897649*2(20,14) 
+0.0213924*2(22,14) -0.0039392*2(24,14) +0.0006007*2(26,14) -0.0000717*2(28,14) +0.0000090*2(30,14) 

0.0001994 : 3: +0.0387722*2(3,3) +0.1746705*2(5,3) +0.3544766*2(7,3) +0.2795820*2(9,3) -0.2630382*2(11,3) 
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-0.4046044*Z(13,3) +0 . 6071129*Z(15 ,3) -0.3742479*Z(17,3) +0. 1447577*Z(19,3) -0 . 0399704*Z(21 ,3) 
+0.0084474*Z(23,3) -0.0014183*Z(25,3) +0.0001974*Z(27,3) -0.0000221*Z(29,3) +0 . 0000026*Z(31 ,3) 

0.0001963 : 1: +0 . 0039977*Z(1 , 1) +0 . 0501564*Z(3, 1) +0. 1886747*Z(5, 1) +0 . 3560673*Z(7, 1) +0 . 2581983*Z(9 , 1) 
-0.2814148*Z(11,1) -0.3836754*Z(13,1) +0. 6083843*Z(15 , 1) -0.3837324*Z(17, 1) +0 . 1508905*Z(19 , 1) 
-0.0422465*Z(21,1) +0 . 0090396*Z(23, 1) -0 . 0015357*Z(25 , 1) +0.0002159*Z(27, 1) -0 . 0000245*Z(29 , 1) 
+0.0000029*Z(31,1) 

0.0001848 : 8: +0.2521278*Z(8,8) +0.4396161*Z(10,8) -0.0455476*Z(12,8) -0.5579080*Z(14,8) 
+0.5728641*Z(16,8) -0.2999791*Z(18,8) +0. 1045395*Z(20,8) -0.0267374*Z(22,8) +0.0053315*Z(24,8) 
-0.0008530*Z(26,8) +0.0001150*Z(28,8) -0 . 0000121*Z(30 ,8) +0.0000016*Z(32,8) 

0.0001789 : 11: -0 .4461735*Z(11 , 11) -0.2838305*Z(13, 11) +0.6602212*Z(15, 11) -0.4854523*Z(17, 11) 
+0.2108168*Z(19,11) -0.0643653*Z(21,11) +0 . 0148112*Z(23 , 11) -0.0027123*Z(25, 11) +0 . 0004011*Z(27 , 11) 
-0 . 0000515*Z (29 ,11) +0 . 0000047*Z (31 ,11) 

0.0001630 : 15: +0 . 7317763*Z(15, 15) -0.6045873*Z(17, 15) +0 . 2977377*Z(19 , 15) -0 . 0984783*Z(21 , 15) 
+0.0246634*Z(23,15) -0 . 0048055*Z(25 , 15) +0 . 0007755*Z(27 , 15) -0.0000998*Z(29, 15) +0 . 0000127*Z(31 , 15) 

0.0001621 : 6: +0. 1321892*Z(6,6) +0.3482620*Z(8,6) +0 . 3088642*Z (10,6) -0.2477773*Z(12,6) -0.4081762*Z(14,6) 
+0.6063296*Z(16,6) -0.3809198*Z(18,6) +0 . 1519987*Z(20 ,6) -0.0436556*Z(22,6) +0.0096556*Z(24,6) 
-0.0017060*Z(26,6) +0 . 0002504*Z(28,6) -0.0000300*Z(30,6) +0.0000036*Z(32,6) 

0.0001501 : 4: +0.0559858*Z(4,4) +0.2088990*Z(6,4) +0 . 3571657*Z(8 ,4) +0. 1959894*Z(10,4) -0 . 3351164*Z(12 ,4) 
-0.2950686*Z(14,4) +0.6047083*Z(16,4) -0.4289290*Z(18,4) +0. 1863726*Z(20,4) -0.0574666*Z(22,4) 
+0.0135356*Z(24,4) -0 . 0025371*Z(26 ,4) +0 . 0003927*Z(28 ,4) -0.0000500*2(30,4) +0.0000060*Z(32,4) 
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